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Abstract 

The decision-making trial and evaluation laboratory (DEMATEL) method employs 
expert assessments expressed by crisp values to construct a group initial direct-rela- 
tion (IDR) matrix. However, it tends to be a low-precision expression, especially in 
complex practical problems. Although significant efforts have been made to improve 
the DEMATEL method, these improvements tend to neglect individual characteris- 
tics and group consensus, resulting in unconvincing decision results. This study pro- 
vides a Dezert-Smarandache theory-based group DEMATEL method with reach- 
ing consensus. In order to reasonably determine the group IDR matrix, basic belief 
assignment function is employed to extract expert assessments and the proportional 
conflict redistribution rule no.5 of DSmT is employed to make fusion to derive the 
temporary group IDR matrix. Moreover, the consensus measures at both expert level 
and pair-factors level are calculated to determine whether the acceptable consen- 
sus level has been reached or not. If the required consensus level is not reached, a 
feedback mechanism will be activated to help experts reach a consensus. A consen- 
sus group IDR matrix for the group DEMATEL can be obtained with the help of 
feedback mechanism, based on which an algorithm is summarized for the proposed 
method to identify major factors in a complex system. Finally, numerical compari- 
son and discussion are introduced to verify the effectiveness and applicability of the 
proposed method and algorithm. 
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1 Introduction 


Between 1972 and 1976, the Science and Human Affairs Program of Battelle 
Memorial Institute of Geneva developed the decision-making trial and evaluation 
laboratory (DEMATEL) method. This method aimed to describe the basic con- 
cept of contextual relations and identify cause-effect chain factors for a complex 
decision problem in an understandable manner by addressing the influence rela- 
tions among factors given by experts (Fontela 1974; Gabus 1973). It was consid- 
ered to be a credible decision-making method. 

The DEMATEL method has been extensively used to solve complex decision 
problems because of its simplicity and effectiveness, including problems per- 
taining to hospital service quality (Shieh et al. 2010), decision making (Michnik 
2013), sustainable supply chain management (Liang et al. 2016; Ren et al. 2013), 
etc. In special, in order to determine the weights of factors by considering their 
relations, the DEMATEL method is also extended into decision making fields, 
such as analytic hierarchy process (AHP) (Sara et al. 2015; Balsara et al. 2019), 
analytic network process (ANP) (Gölcük and Baykasoglu 2016; Lan and Zhong 
2016; Ghaemi Rad et al. 2018; Quezada et al. 2018; Chen et al. 2019; Dinger 
et al. 2019), and technique for order preference by similarity to ideal solution 
(TOPSIS) (Baykasoglu and Gölcük 2017; Baykasoglu et al. 2013). In the DEM- 
ATEL method, the initial decision information is always subjectively given by 
experts in the form of crisp values and calculated to obtain an individual or group 
initial direct-relation (IDR) matrix by simple operations (e.g., weighted sum). 
However, such descriptions and operations are considered to hardly reflect the 
vagueness of the real world (Chen and Hwang 1992). Therefore, scholars have 
carried out some research to improve the DEMATEL combined with fuzzy theory 
(Wu and Lee 2007). Several fuzzy DEMATEL methods have been introduced. 
For examples, Abdullah et al. (2019) introduced interval-valued intuitionistic 
fuzzy numbers to improve the judgement of DEMATEL in a group decision-mak- 
ing (GDM) environment. Addae et al. (2019) used a two-step fuzzy DEMATEL 
method to solve a practical problem. Asan et al. (2018) proposed a new interval- 
valued hesitant fuzzy approach to DEMATEL to explicitly deal with hesitation in 
expert assessments and offered a better representation of uncertainty, etc. (Bhatia 
and Srivastava 2018; Acufia-Carvajal et al. 2019; Cui et al. 2019; Du and Zhou 
2019). 

Obviously, all these extensions have made great contributions to the DEMA- 
TEL method, and its ability of dealing with complex problems can be strength- 
ened to some extent. As currently defined, instead of one expert, a panel gives 
the assessments on the influence relations of factors, and multiple experts arrive 
at an acceptable result (Cheng and Lin 2002). This process makes it necessary to 
extend the traditional DEMATEL method to a group DEMATEL method, which 
belongs to GDM problems. Although some literature considers the DEMATEL 
from the perspective of GDM, we believe that these group DEAMTEL methods 
have a lot of room for further improvement in both expert assessment extraction 
and group IDR matrix construction. 
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Firstly, the expert assessment extraction of group DEMATEL method should be 
improved to obtain accurate decision information in accordance with experts’ cogni- 
tive competence. As introduced earlier, DEMATEL is a GDM method totally based 
on expert assessments to conduct later computation and analyzation. Considering 
the complexity of reality and individual characteristics of experts, there is no doubt 
that not every expert is proficient in all areas. In other words, experts may give 
incomplete or uncertain assessments to the influence relations among factors for a 
specific complex practical problem. However, all of the existing group DEMATEL 
methods default that each expert can give a definite assessment for every pair of 
factors by means of crisp values or fuzzy numbers, and neglect the problem that the 
assessments given by experts may be incomplete in practice. If experts who cannot 
give definite assessments with crisp values or fuzzy numbers are required to give 
the assessments in those forms, the final decision of DEMATEL resulted from the 
ineffective assessments may be erroneous. Therefore, understanding how to depict 
and fuse incomplete information from experts is of great importance to improve the 
properties of the group DEMATEL method. Fortunately, the basic belief assignment 
(BBA) function, as a key concept in the Dezert-Smarandache theory (DSmT) of 
evidence (Smarandache and Dezert 2004), can directly express the uncertainty by 
assigning probability to the subsets of a set composed of multiple objects rather than 
to each individual object (Du et al. 2018c). The BBA functions generated from dif- 
ferent evidence sources (experts) could be well fused by the proportional conflict 
redistribution rule no. 5 (PCR5). All these features exactly meet the requirements of 
the group DEMATEL method. Accordingly, DSmT is used to extract and fuse expert 
assessments to derive the group IDR matrix in this paper. Moreover, the differences 
in experts’ knowledge backgrounds and cognitive abilities on a particular problem 
are reflected by expert weights (Du et al. 2018b). Expert weights are reflected as 
discounting parameters to reflect one’s relative importance in a group during fusion 
process. In this paper, expert weights are calculated based on similarity functions of 
expert assessments. 

Secondly, the group IDR matrix construction of group DEMATEL method 
should be improved to obtain the acceptable decision results in accordance with 
experts’ satisfaction. In the GDM, “group” refers to not only the number of experts 
merely, but also the experts who have common interests in reaching a consensus for 
the ultimate satisfactory results despite individual differences. This principle helps 
to reduce biased evaluations and inherent partiality in GDM processes (Biiytikéz- 
kan and Çifçi 2012). Unfortunately, the group IDR matrix in the existing group 
DEMATEL methods is frequently constructed by making arithmetic average values 
for individual IDR matrices, while whether experts agree with the group results is 
scarcely concerned. The group IDR matrix plays a fundamental role in the entire 
DEMATEL processes and it has a significant influence on the effectiveness of final 
results. If strong inconsistency and conflict exist among experts, the group IDR 
matrix may not be able to precisely describe the real influence relations of factors. 
Therefore, it is particularly important to construct a group IDR matrix according to 
the consensus rules, that is, to construct a group consensus IDR matrix by reaching 
general or widespread agreement among the experts involved in the GDM processes 
(Pérez et al. 2018). Fortunately, Herrera-Viedma proposed a rational consensus 
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model in GDM composed of a selection process and a consensus-reaching pro- 
cess (CRP), which had become a hot issue in the recent GDM area. The CRP has 
been successfully introduced to make GDM with different situations, such as hesi- 
tant fuzzy preference relations, Delphi processes, multi-attribute large-scale GDM, 
sentiment analysis, and virtual reality industry (F. Herrera et al. 1997; Liang et al. 
2016; Michnik 2013; Ren et al. 2013; Shieh et al. 2010). Traditionally, unanimous 
agreement of all experts in CRP is required. However, the desired result can hardly 
be achieved because of the diversity of opinions, knowledge, and experiences of 
experts. Therefore, the concept of “soft consensus” has been provided, in which, 
“soft” means better reflecting all possible levels of agreement by setting an accept- 
able consensus level (CL) threshold value (such as 0.8 rather than 1) and guiding the 
CRP until a high-level agreement is achieved among the individuals. Soft consen- 
sus can be reached through an iterative dynamic process with several collection and 
adjustment rounds (Kacprzyk and Zadrozny 2010; Cabrerizo et al. 2017). Hence, in 
this paper, we consider a soft CRP in the group IDR matrix construction, and trans- 
form the original static group DEMATEL problems into dynamic ones. 

The motivation of this paper is to improve group DEMATEL method according 
to the following three aspects. Firstly, the BBA function is used to extract expert 
assessments to accurately express uncertainty and incompleteness, thus reducing 
the loss of accuracy. Secondly, the initial assessments are discounted with expert 
weights by using Shafer’s discounting method. Moreover, the PCRS of DSmT is 
used to fuse the discounted assessments to overcome the defects in the intuitional 
paradox of Dempster’s combination rule. Thirdly, we apply a soft CRP to help reach 
an acceptable CL in the construction of group IDR matrix to ensure the consistency 
and satisfaction among experts. This paper is organized as follows. Section 2 briefly 
introduces the basic knowledge of DEMATEL and DSmT. In Sect. 3, the DSmT- 
based group DEMATEL method with reaching consensus is proposed, and the cor- 
responding algorithm is summarized. In Sect. 4, the numerical comparison and dis- 
cussion are provided to demonstrate the performances of the proposed method and 
algorithm. In Sect. 5, the conclusion is drawn and the future research directions are 
briefly discussed. 


2 Preliminaries 


In order to facilitate the later formulation, some basic concepts of DEMATEL and 
the Dezert-Smarandache theory (DSmT) are given in this section. 


2.1 DEMATEL Method 


The DEMATEL method is an effective way to analyze the influence relations among 
factors of a system. Through an analysis of the total influence relation of factors by 
the DEMATEL method, we can obtain a better understanding of structural relations 
and an ideal way to solve complicated system problems. Consider that a group of 
experts are invited to assess the influence relations for a set of factors F = {f,,...,f;} 
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with a set of grade levels {0, 1, 2}, where the expert set is E = {e}, ... , eg }, and the 
degree of influence to which he or she believes factor f; has an effect on factor f; 
(denoted by f; > f;) is assessed by expert e, and denoted as 8i Vi, j,k. The meaning 
of each element in the set of grade levels is that O denotes “no influence”, 1 denotes 
a “low influence” and 2 denotes a “very high influence”. The central concepts in the 
DEMATEL are defined as follows. 


Definition 1 (Singhal et al. 2018) Suppose a pairwise comparison of influence 
degree from the ith to the jth factor given by the kth expert e, is denoted as Bi, with 
0-2 grade levels, and the grade levels given by each expert form a L x L non-nega- 
tive answer matrix G* = [si] xp K = 1,...,K. The group IDR matrix, which repre- 
sents the initial direct relation between each pair of factors derived from experts, is 
obtained by calculating the average values of all experts’ answer matrices as 
G = [gyl where gy = Di, 8 /Kij = 1,...,L. 

In Definition 1, the 0, 1, 2 grade levels mean “no influence,” “medium influ- 
ence,” and “high influence,” respectively. Note that, the diagonal elements of each 
answer matrix G* are all set to zero, which means that the factors do not influence 
themselves. 


99 66 


Definition 2 (Singhal et al. 2018) Suppose the maximal row-wise and column-wise 
sum of matrix G is g’ = max{ max) <<, De Zij MAX <j<L E sy): then the nor- 


malized IDR matrix D = [dij] xz Can be computed according to Eq. (1). 


D=G/ g (1) 


Definition 3 (Singhal et al. 2018) Suppose the direct and indirect relations among 
several factors are represented by the total relation matrix, and it is defined as in 
Eq. (2). 


A= lim (D +D? +- +D") = DU-Dy! 2) 


Some kinds of extensions are further discussed to strengthen the original DEMA- 
TEL. One kind of extensions is used to overcome the drawback that raising the nor- 
malized IDR matrix to the power of infinity may not converge to zero, and hence, 
the total relation matrix may not converge (see Eq. (2)). A very small positive num- 
ber u (e.g., u = 1075) is introduced in the maximal row-wise and column-wise sum 
of matrix G as g” = max (max, sir 2i Zi H + MAX) <j<L EL sy}: Other steps 
remained unchanged as in the original DEMATEL. The revised DEMATEL guaran- 
tees that the normalized IDR matrix to infinite power will converge to zero and that 
the total relation matrix can be obtained smoothly (Lee et al. 2013). 


Definition 4 (Singhal et al. 2018) Suppose r and c represent the sum of rows and the 


sum of columns of the total relation matrix A. According to A = [aj]; r and c can 
be defined as follows: 
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(3) 


where r; shows the total influence, both direct and indirect, given by the factor f; to 
other factors; c; shows the total influence, both direct and indirect, received by the 
factor f; from other factors; r; + c; is defined as the prominence, showing the degree 
of the important role that the factor f; plays in the complex system; and r; — c; shows 
the net influence that the factor f; contributes to the complex system. Note that, if 
r; — c; is positive, the factor f; is a net causer; if r; — c; is negative, the factor f; is a 
net receiver. 


2.2 Dezert-Smarandache Theory 


DSmT, jointly proposed by Smarandache and Dezert (2006), can be used to obtain 
more accurate fusion results of BBA functions especially in high conflicting infor- 
mation cases. It has a series of proportional conflict redistribution rules to make 
fusion for evidences (Du et al. 2019), among which, PCR5 is the most widely used 
one with the advantages in dealing with conflict belief functions. For example, it 
provides the appropriate redistribution of conflict beliefs and can produce a reason- 
able fusion result even in highly conflicting cases. These attractive features motivate 
the use of DSmT in GDM problems, such as map reconstruction of robot (Singh 
et al. 2008), decision making support (Huang et al. 2014; Liu et al. 2011), target 
type tracking (Faux and Luthon 2012), image processing (Liu et al. 2012a), data 
classification (Lian et al. 2015; Liu et al. 2013, 2015a, b, 2016), clustering (Denceux 
et al. 2015; Liu et al. 2012b, 2015b), and so on. 

In DSmT framework, the frame O = {0}, ...,0y} is a finite set of Y exhaustive 
propositions that are not necessarily mutually exclusive. The hyper-power set D® is 
defined as the set of all composite propositions built from elements of © with U and 
N operators, such that (Guo et al. 2016): 

(i) Ø, 0),...,0y € D®; 

(ii) If O,, 0y E D®, then 0, U 0y E€ D® and 0, N Ay € DO; 

(iii) No other elements belong to D®, except those obtained by using rules (i) or 
(il). 

Definition 5 (Smarandache and Dezert 2006) Suppose © = {6,,...,0)} is a set of 


exhaustive propositions; then the basic belief assignment is defined over the hyper- 
power set D®. If the mapping function m : D? = [0, 1] could fulfill the following: 
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mø) = 0, DY m(9) = 1 (4) 
OCD? 
then m(-) is called the BBA function. If m(@) > 0, 0 is called a focal element. 
In DSmT framework, PCR5 for making fusion for two pieces of evidence is intro- 
duced as follows. 


Definition 6 (Smarandache and Dezert 2006) Suppose the BBA functions of two 
pieces of evidence are m, and m, on D®; then, PCR5 to fuse m; and m, can be 
defined as follows: 





6 2 el” 6 2. o" 

Y momo J [= AATAL 
onoo onapo L™m(0)+m,(0™)  m(0) + m, (0) é 

mpcrs(9) = 4 gr gnep ono" =g , OED" andd FD 


0, 0=9 


(5) 
3 The Proposed Method 


In the proposed method, we construct the group IDR matrix in the DEMATEL by 
following two steps. Firstly, we develop an expert assessments extraction and fusion 
mechanism with DSmT to obtain a temporary group IDR matrix (see Sect. 3.1). Sec- 
ondly, we activate the soft CRP based on the current individual and group IDR matri- 
ces to reach a soft consensus among experts (see Sect. 3.2). 


3.1 Expert Assessment Extraction and Fusion 


In traditional DEMATEL, experts are advised to give their assessments of the influence 
relations among factors by means of crisp values to construct the group IDR matrix 
G= [gil ix (see Definition 1), which is a rough extraction method with low preci- 
sion. When restricted by expert experiences and knowledge, the assessments given by 
experts may be ignorant and partially credible. Taking the influence degree for f; > f; 
and expert e, as an example, the expert may well know this problem and he or she 
can give a definite result by using one of the given grade levels {0, 1, 2}. On the con- 
trary, the expert may not know this problem at all, in which case he or she cannot give 
any assessment information about the influence degree. In most instances, the expert 
knows this problem to a certain extent and he or she may point out that the influence 
degree belongs to several grade levels but may not be sure which is the best one. Obvi- 
ously, experts cannot give their assessments by means of crisp values in the latter two 
situations—that is, the original extraction method is not practical and does not consider 
experts’ personality. Thus, we employed the BBA function to extract expert assess- 
ments in this paper as shown in Fig. 1. 

As shown in Fig. 1, the expert assessment extraction and fusion mechanism includes 
three major steps. (1) Extracting expert assessments with BBA function; (2) Discount- 
ing these assessments with Shafer’s discounting method; (3) Fusing the discounted 
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assessments with PCRS. Next, we will give a detailed definition and explanation of the 
procedures involved in the above processes. 

Let the frame of discernment be © = {0,, 05,03} = {0, 1,2}, which can be seen as 
the discernment frame of DSmT. Expert e, is asked to assess the influence degree for 
Ji > f; and his/her assessment is allowed to be expressed by the BBA function as given 
in Definition 5 and shown in Eq. (6). 


b; = fu b0) È, By(8) = 1; Bi(8) > 0, 0 C 0, vija) (6) 
oco 


All of the BBA functions make up the individual IDR matrix for expert ep, 
denoted as B* = [bile with i,j = 1,...,L and k= 1,...,K. Note that the basic 
beliefs in bi can be assigned not only to singleton grade levels but also to any sub- 
sets of ©, thereby it is allowed such an assessment (also called a piece of evidence) 
to be profiled by a BBA defined on the hyper-power set D®. It is capable of reflect- 
ing ignorance in expert assessments, and the basic beliefs in bi can be given to © 
(global ignorance) or to 8 C © (local ignorance) according to the unknown and par- 
tial assessments (Du and Xu 2017). Expert assessments could be expressed precisely 
through BBA functions as discussed, laying a great foundation for later DEMATEL 
procedures. 


Example 1 Assume expert e} points out the influence degree for f, > f; has 30% 
probabilities belonging to 0, and 70% probabilities belonging to 0,. Thus, his or her 
assessment is described as bi() ={(9,, 0.30), (02, 0.70)}. Expert e, points out the 
influence degree has 20% probabilities belonging to 0, and has 80% probabilities 
belonging to (0, U 03) but is not sure which is the best one. Therefore, his or her 
assessment is described as b;(8) = {(0,, 0.20), ((@ U 83), 0.80)}. Expert e, points 
out the influence degree for f; > f; has 100% probabilities belonging to (0, N 03). 
Thus, his or her assessment is described as b30) = {((8. N 03), 1.00)}. Expert e, 
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cannot give any information about the influence degree. Therefore, his or her assess- 
ment is described as b30) = { (©, 1.00)}. 


Normally, expert weight is subjective and relative to reflect the importance of 
one’s assessments in a group, and it is usually denoted by w in [0, 1], with O and 
1 respectively standing for not important at all and the most important (Du and 
Wang 2017; Du et al. 2018a). It is obvious to find that the expert weight can 
be determined by AHP (Saaty 2003), ANP (Saaty 2007) or Delphi (Stebler et al. 
2015), and it can also be determined subjectively according to the requirements 
of actual issues (Li et al. 2015). As discussed, expert weight is used to account 
for one’s relative importance among all experts—that is, the closer one’s assess- 
ments are to others’ assessments, the more important the expert is likely to be 
(Huang et al. 2014). Hence, in our opinion, expert weight could be calculated 
indirectly based on the similarity between one expert and other experts—that is, 
the higher the similarity of assessments between the expert and others, the larger 
the weight of the expert. Expert weight is directly proportional to the similar- 
ity between expert assessments and others’ assessments. Thus, we use evidential 
distance to depict the similarity of expert assessments, based on which we could 
calculate expert weight in simple ways. The Euclidean evidential distance and 
Euclidean evidential similarity, which have little computation complexity and fast 
convergence speed, are defined as follows. 


Definition 7 (Li et al. 2011) Suppose m, and m, are two BBA functions for the same 
frame of discernment O, 0, is the n th element of D®, and |D°| is the cardinality of 
D®. The distance between m, and m, is defined as follows. 


Distg(m,, m) = L yr [m,(0,,) — m (0,)]? (7) 


V2 n=1 


Definition 8 (Li et al. 2011) Suppose m and m, are two BBA functions for the same 
frame of discernment O, 0, is the n th element of D®, and |D°| is the cardinality of 
D®. The Euclidean similarity function Sim,(m,,m,) is defined based on the Euclid- 
ean evidential distance as follows. 


. 1 g i 
Simp(m,,m,) = 1 — R po [m,(6,,) — m,(6,,)]? (8) 


Example 2 Assume two pieces of evidence, b, ={(6,,1)} and 
b, = {(0,, 0.4), (@,,0.6)}, and that the distance between them could be computed 
by Eq. (7): Dist,(b,,b,) = y0 -0.4)2 + (1 - 0.6)2//2 = 0.40. The similarity 
between them is computed by Eq. (8): Sim,(b,, b2) = 1 — 0.40 = 0.60. 

On the basis of the previous definitions, the similarity between experts can be 
directly calculated. Thus, we define the following expert weight calculation method. 
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Definition 9 (Yang and Wang 2010) Suppose the individual IDR matrix consisting 
of all influence relations among factors given by expert e, is B*= [biex with 
ij=1,...,L and k= 1,...,K. The similarity between experts e, and e, on each 
pair of factors f; > fj is Sim,(b, bk ), and the similarity between any two experts is 
computed as s = peer 4 Simp (bi, bt) /L?, where L? denotes the quantity of fac- 
tor pairs, making up the similarity matrix as follows. 


1 glk na. glk 
S= [lace =| soe 1 oh (9) 
sK! kk 1 


The support of e, can be obtained by adding all of the elements in the simi- 
larity matrix S that are related to expert e, except for self-similarity, i.e., 
Sup(e,) = seer sk =1,...,K, where Sup(e,) represents the support degree 
of expert e, received from other experts. By normalizing them, we could obtain the 
expert credibility Crd(e,) which is generally regarded as expert weight w* as follows. 


K 


wk = Crd(e,) = Sup(e,)/ È, Sup(e,), k = 1,...,K (10) 
k=1 


Example 3 Assume experts ¢€,,€,€3 give their assessments for the influence rela- 
tions among factor set F = {f,,/5,f;}, respectively, as follows: 


{(0,, 1.0)} {(03,1.0)} — {(83, 1.0)} 
B! = {(@, 0.7), (83, 0.3)} {(0, 1.0)} {(0,, 0.4), (3, 0.6) } 

{(0;, 1.0)} {(02,1.0)}  {0,, 1.0)} 

{(,, 1.0)} {(0, 1.0)} {(03, 1.0)} 
B? =| {(6,,0.8),(0,,0.2)}  {(04, 1.0)} {(0,, 1.0)} 

{(0,, 1.0)} {(0;,0.5),(8),0.5)}  {(0,, 1.0)} 

{(0,, 1.0)} {(0,1.0)}  {(01, 0.2), (0; U 03, 0.8)} 
B? =| {(0,U63,1.0)} {(0,1.0)}  {(84, 1.0)} 

{(0,, 1.0)} {(0;,1.0)}  {(0,, 1.0)} 


The similarity matrix can be calculated by Definition 9 as follows: 
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1.00 047 0.29 
S=| 0.47 1.00 0.41 
0.29 041 1.00 


We calculate the support degree of experts as Sup(e,) = 0.47 + 0.29 = 0.76, 
Sup(e,) = 0.47 + 0.41 = 0.88, Sup(e3) = 0.29 + 0.41 = 0.70. Then we compute the 
expert weight as w! = Crd(e,) = 0.76/(0.76 + 0.88 + 0.70) x 
{01 Pi) (02, P3)» (03, PD} w? =0.87/ 2.34 = 0.37, w° = 0.70/ 2.34 & 0.30. 

As mentioned above, expert weight should be considered through Shafer’s dis- 
counting method, which multiplies the masses of focal elements by the expert 
weight and transfers all of the remaining discounted mass to the full ignorance ©. 
Mathematically, Shafer’s discounting method can be given as in Definition 10. 


Definition 10 (Shafer 1996) Suppose the BBA function is b as in Eq. (4), and w is a 
parameter to discount the evidence, 0 < w < 1. Then, Shafer’s discounting method is 
defined as follows: 


0, 0=ğ 
m0)=<+w-b@),0dc 0, 046 (11) 
w:-b(@)+l—w,9=0 


If the sum of weights is equal to 1, then the discounting parameter is usually 
derived by standardizing the weights as w* = w‘/max(w*|k = 1,...,K). The dis- 
counted expert assessments could be described as follows. 


m; = fu m0) 2 = 1;m;(0) = 0,0 C ovije} (12) 
ac 


The discounted individual IDR matrix is denoted as M* = [mi] ix, With 
ij=1,...,Landk=1,...,K. 


Example4 Assumeexpert weightisW = {w! = 0.40, w? = 0.20, w° = 0.20, wt = 0.20}. 
The assessments given by four experts are the same as in Example 1. The discounting 
parameters for the four experts can be standardized as w= w!/ max(w!, nts yw )= 
0.40/0.40=1.00, and similarly we get W = W° = W° = 0.20/0.40 = 0.50. Taking them 
into Eq. (11), the discounted assessments are as follows: m;,() = {(0,, 0.30), (03, 0.70)} 
(8) = {(0,, 0.10), (0 U 43), 0.40), (O, 0.50) }, m;(8) = {(8) N 03, 0.50), (©, 0.50)}, and 
m9) = {(O, 1.00)}. 


The DSmT framework with PCR5 rule can be used to make fusion for individ- 
ual discounted assessments as in Eq. (12) and the fusion result can be described as 
follows: 
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oco 
The determined group IDR matrix is denoted as M = [mj]7y7. 


Example 5 Assume two pieces of evidence are: 
m;(8) = {(0,, 0.20), (82, 0.30), (O, 0.50) }m (0) = {(0,,0.70), (O,0.30)}. Taking 
them into Eq. (5), the fusion results are obtained as m,(@) = {(84, 0.09), (82, 0.76), 
(0, 0.15)}. 


3.2 The Soft CRP for Group DEMATEL 


Soft CRP can be designed with the aim of supporting experts until a group con- 
sensus is reached by following several discussion and adjustment rounds. As dis- 
cussed in Sect. 1, situations in which all of the experts agree with each other 
unanimously are rare or not desirable in the decision making process. Since 
the “unanimous consensus” of conflict tolerance, which has the ability to sat- 
isfy all pairs of BBA functions, rarely exists, the choice of a “soft consensus” 
is largely subjective and application oriented (Liu 2006). That is, setting an 
acceptable CL threshold value (suppose that the value is 0.8 there) to guide the 
whole process to reach group consensus. Basically, this process relies on mak- 
ing assumptions about experts’ willingness to change their opinion or preferences 
(Herrera-Viedma et al. 2007; Pérez et al. 2011). Initially, consensus measures are 
computed based on individual and group IDR matrices to determine whether an 
acceptable CL has been reached or not. If so (CL>0.8), the soft CRP is finished 
and a consensus group IDR matrix can be obtained. Otherwise (CL <0.8), the 
feedback mechanism will be activated, and the experts who are not contributing 
to the consensus are identified and the advice about how to alter their assessments 
is generated (Xu and Wu 2011). 

This consensus measure can indicate the current consensus situation through- 
out the soft CRP. According to the characteristics of group DEMATEL, the con- 
sensus measure is categorized into two levels, i.e., pair-factors level and expert 
level. The pair-factors level is the most basic level and reflects the original con- 
flict degree. At this level, experts make some adjustments to increase consistency 
among the group. The expert level is the highest level, which reflects the conflict 
degree between a specific expert and the group as a whole for the collected results 
on all pairs of factors (Wu et al. 2017). We employ these two levels of consensus 
measures to identify inconsistencies between experts and group. 

The calculation methods for consensus measures can be derived based on the 
similarities between these assessments. At the pair-factors level, the consensus 
measure values can be calculated based on the Euclidean similarity function as in 
Definition 11. At the expert level, the consensus measure values can be calculated 
by adding the consensus measure values on all of pair-factors for an expert as in 
Definition 12. Obviously, the order for calculating the two levels of consensus 
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measure is pair-factors level at first and then expert level. When determining the 
inconsistency assessments that need to be modified, the order is the highest level 
at first and then the most basic level. Experts only adjust the conflicting assess- 
ments based on the recommendations at the pair-factors level. 


Definition 11 In the ¢ th round, suppose the Ass sormieni of the pair of factors f; > f; 
given/calculated by expert e, and group are be and mi. Then, the consensus meas- 


ure values can be calculated by the Euclidean similatity function as follows. 


a4 av au | = mi,(,)P (14) 


Definition 12 In the ż th round, suppose the consensus measure value of expert e, 
for f; > f; at pair-factors level is e t Vi, j,k, then the consensus measure value of 


expert e, at expert level in this fond can be calculated as follows. 


L 


cht = > ee (15) 


ij= Lit 
where L? denotes the quantity of factor pairs. 


The closer the value c is to 0, the greater the conflict degree; the closer the value is 
to 1, the smaller the conflict. Obviously, c“’ = 0 indicates complete conflict between 
expert e, and group, and c* = 1 indicates no conflict between them. If each expert’s 
consensus measure c*" is larger than the acceptable CL threshold value, indicating that 
existing conflict degree can be accepted by all experts and group consensus is reached, 
then the soft CRP is finished and the current collected group IDR matrix M‘ = [mi] LXL 
is taken as a consensus result. Otherwise, if one’s consensus measure c* is lower than 
the acceptable CL threshold value, indicating that strong conflict degree among experts, 
then the feedback mechanism is carried out to help the inconsistent experts adjust their 
assessments to enhance group consensus, and a new round (t + 1) of group IDR matrix 
construction should be initiated. 

Note that, expert weights may be not fixed during several of the rounds for the rea- 
son that expert assessments may be changed in the context of not reaching the accept- 
able CL. Thus, in each new round, we need to recalculate expert weights based on the 
latest assessments. In order to prevent the collective assessments from failing to con- 
verge after several discussion rounds, we incorporate a maximum number of rounds 
(tyax) in the soft CRP to develop. It can be ensured that the feedback mechanism will 
not be carried out when t = ty,x and the current collected group IDR matrix will be 
taken as the final result even if the acceptable CL has not been reached yet. 

Suppose the acceptable CL threshold values at expert level and pair-factors level are 
£ and £p. Then the processes of feedback mechanism can be divided into the following 
two steps. 
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Step 1 Experts whose consensus measure values at expert level are lower than the 
threshold value £, in the ¢ th round are identified as follows. 


È, = {k| <e,} (16) 

Step 2 For the identified experts in step 1, their assessments for such pairs of factors 

that the consensus measure values are lower than the threshold value Ey are identified 
as follows. 


F, = {k,i > jlc% < ep ak € È} 7) 


When E and F have been identified, personal advices for experts will be generated 
to reduce the conflict caused by the moderator. Since pair-factors level is the most basic 
level and contains the original assessments of experts, only the experts whose consen- 
sus measure values c* are lower than the threshold value £, may obtain the adjusted 
advices. After the identified experts all finish their adjustments, a new round of group 
IDR matrix construction will be carried out to determine the consensus situation. It is 
obvious to find that the above procedure is very similar to the Delphi method. 


3.3 DSmT-Based Group DEMATEL Method and Algorithm 


The construction of a group IDR matrix includes two main processes. The first is 
collection process, which focuses on expert assessment extraction and fusion; and 
the second is soft consensus reaching process, which aims at reaching an acceptable 
CL among experts. After the above two processes, the final group IDR matrix M can 
be constructed. Obviously, M is different from the original IDR matrix G as given 
in Definition | for two reasons. (1) M is a consensus group IDR matrix that satisfies 
all experts, while G does not consider consensus and satisfaction of experts. (2) M 
is made up of BBA functions that is capable of reflecting local or global ignorance 
in experts’ mind, while G is made up of crisp values. Note that, the IDR matrix 
included in the DEMATEL is required to be exact influence degrees rather than 
BBA functions. Therefore, we apply the generalized pignistic probability as in Defi- 
nition 13, which is an extension of the original pignistic probability (Smets 2005), 
to reassign the local and global ignorance in the BBA functions to singleton grade 
levels. Then, we calculate the expected value to derive the exact influence degrees. 


Definition 13 (Smarandache and Dezert 2009) Suppose the frame of discernment in 
the DSmT framework is O = {6,,..., Oy}, then the generalized pignistic probability 
for VA € D® can be calculated by Eq. (18). 





A 
pras Y EA no us 


where |X| denotes the cardinal of proposition X. 
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Fig. 2 The process of proposed method 


Taking m, in the final group IDR matrix M into Eq. (18), the transformation 
results can be obtained and recorded as {(6,, P}) ...,(Oy, P;,)}, where P} denotes the 
probability of influence degree 0,,y = 1, ... , Y. The subscript Y is equal to 3 for the 
reason that the frame of discernment is defined as O = {01,02,03} = {0, 1,2} in this 
paper. Which influence grade level or influence degree is attached to the relationship 
of fi > f;? We follow two kinds of principles to solve this problem: (1) Highest 
probability principle. The grade level with the highest probability 2; = 0, is chosen 
as the final influence grade level, where {9,|P; = max(P;., ws Pb (2) Expected 
value principle. The expected value g, = 2a 0, xP; Vi, j is calculated as the final 
influence degree. According to the real situation, we obtain the final influence 
degree for fi > f; by using one of the two principles, and then the final group IDR 
matrix G = [8;]zxz can be constructed. 

The process of DSmT-based group DEMATEL is illustrated as in Fig. 2. As 
shown in Fig. 2, the complete steps of DSmT-based group DEMATEL can be sum- 
marized as follows. 

Step 1 Define group DEMATEL problem and soft CRP parameters Suppose the 
set of factors is F = {f1, ....f,}, the set of experts is E = {e),...,ex}, the set of 
grade levels is © = {6,, 05, 03}, and the threshold value to filter out major factors is 
n. Then, the acceptable CL threshold values at expert level and pair-factors level are 
£, and £p, the round counter is f (t is set to one at first), and the maximum number 
of rounds is ty,ay. A moderator is invited to participate in the group decision making 
process and responsible for managing the whole process such as consensus measure 
value calculation, inconsistent expert identification, the maximum number of round 
judgement. 
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Step 2 Expert assessments extraction and fusion Experts are advised to indicate 
the influence degree to which he or she believes factor f; has an effect on factor f; 
(denoted by f; > H) under the framework of DSmT. The BBA functions assessed for 
each pair of factors f; —> f; by expert e, are described as bi, (as in Eq. (6)), and make 
up an individual IDR matrix B£ = [bile On the basis of expert assessments, 
expert weights W = {w!,...,w} are calculated by Definition 8. Initial assessments 
b; are discounted by Shafer’s discounting method with expert weight as in Defini- 
tion 10 to obtain the discounted assessments m Then PCRS is applied to fuse those 
BBA functions to obtain group assessments m 
M= [mj] LXL 

Step 3 Soft consensus reaching process The moderator calculates the consensus 
measure values oF and c* by means of expert and group IDR matrices to identify 


j Making up group IDR matrix 


whether an acceptable CL is reached in the current round. If so, the soft CRP is fin- 
ished and a consensus group IDR matrix M has been obtained. Generalized pignistic 
probability is introduced to deal with global ignorance and local ignorance as in 
Eq. (18) and the influence degree for f; > f; (Wij) is determined by one of the two 
principles as mentioned above to construct G = [8,];,;, and proceeding to step 4. 
Otherwise, let t = t + 1 and carry out the feedback mechanism. The moderator iden- 
tifies the experts who differed strongly from the group by Eq. (16) and provides 
advices for them by Eq. (17). Then, return to step 2 and follow the same two pro- 
cesses. Note that the moderator also should detect whether the maximum number of 
rounds (means by t= ty,y) has been reached before carrying out the feedback 
mechanism. If so, the mechanism will be ceased, taking the current group IDR 
matrix M as the final result and applying the generalized pignistic probability to 
obtain G = [&; lrx; Then proceed to step 4. 

Step 4 Calculate the normalized IDR matrix To guarantee that the normalized 
IDR matrix to infinite power will converge to zero and that the total relation matrix 
can be smoothly obtained, 8, =max(max, <i<z EL Žij H + MAX jez, ) is 
introduced as shown in Eq. (1) to calculate the normalized IDR matrix D = [d;] LXE 
where d;; = Byli Vi, j. 


L 2 
j=1 ij 


Step 5 Compute total relation matrix The total relation matrix is calculated by 
A=D((I- DY! = [ay ]rx- The influencing degree of f, is computed by r; = Xa aij, 
the influenced degree of f; is computed by c; = pee a,;, the prominence degree of f; 
is computed by r; + c;, and the net influence degree of f; is computed by r; — c; 

Step 6 Obtain the major factors with the threshold value To simplify the com- 
plexity of a system to a manageable level, negligible factors should be filtered out. 
Only those factors whose prominence degrees are greater than the threshold value 
(r; + c; = y) should be chosen and considered in the complex system. Generally, the 
cause-effect relation diagram can be made for the major factors based on their prom- 
inence degrees and net influence degrees. 
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Through the analysis of total relation of the factors by DEMATEL, a better under- 
standing of the structural relation and an ideal way to solve complicate system prob- 
lems can be obtained. The algorithm can be summarized to help moderator manage 
the whole processes. 





Algorithm 1 The algorithm for DSmT-based group DEMATEL method with reaching consensus 
Inputs: Set of factors F={f,...f,}, set of grade levels O={6,,0,,0,} , set of experts E={e,,...,e¢,}, DEMATEL threshold value 77, 
threshold values at expert level and pair-factors level £, and £, , round counter t, accepted maximum number of rounds tyay . 
Outputs: Set of major factors F={f,, ...,f,}- 
Begin 
fl, F={k ir j|Yk,Yi Yj}, by ={G,)} for Vi, j,k 
Step 1: Collection process 
For i=1 to L 
For j=1 to L 
If i#j 
Fork=-ltoK | 
If {k,i—> je F Then 
Give/modify the influence degree bj” by expert e, 
Else ea 
ie le 
Endif 
EndFor 
Calculate expert weight W, ={w"" |k=1,..,K} by Eqs.(8)-(10) 
Calculate mj" in Eq.(12) by discounting b“ with w“ 
Make fusion to get m, =m; ®...®m;" with Eq.(5) 
Else 
Let mi ={(6,,1)} 
EndIf 
EndFor 
EndFor 
Get temporary group IDR matrix in the sth round M‘=[m)],,, 





Step 2: Soft consensus reaching process 
For i=1 to L 
For j=1 to L 
Calculate consensus measure value at pair-factors level by c$“ =1-4/) 
EndFor 
EndFor 
Calculate consensus measure value at expert level by c* =, cj" [È 
Derive inconsistent expert set by Eq.(16) and derive £, = {k| c” < é} 7 
Derive inconsistent pair of factors set by Eq.(17) and derive F = {k,i j|ch <£, Ake E} 
If +Ø and tst,,, Then 
For each {k,i> j}e É, 
Give the advices that expert e, is suggested to make corrections for bj" 
EndFor 
Let @t+1 
Return to Step 1 
Else 
Get final BBA-formed group IDR matrix M‘ =[m;],,, 
Turn to Step 3 
EndIf 


[o*l 


[b (8,) - m;(8,)P ha 


Step 3: Group DEMATEL process 
For i=1 to L 
For j=1 to L 
Take m, into Eq.(18) and get its generalized pignistic probability {(6,, P;), (8,, P7), (0, P?)} 
Use the highest probability or the expected value principle to obtain the final influence degree g, 
EndFor 
EndFor s 
Construct final group IDR matrix GĜ=[8;]na L y 
Calculate normalized IDR matrix by D =G / max(max Ža Š; 4+ max 2 485) 
Calculate total-relation matrix by „A= (a;)n4 = Dir by one. 
Initialize set of major factors by F =Ø 
For i=1 to L 
Calculate the important degree of f, by ,+¢,=))',4,+> a; 
If r+¢,27 se hag he a 
Add f, into the set of major factors by F=FUf, 
EndIf 
EndFor 
End 
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Table 1 Assessments of expert e, in the first round 
hoh Ah Ah irk hoh hoh hoh fo >h fh 


0, 0.00 0.50 0.00 0.00 0.00 0.00 0.50 0.00 0.00 0.00 


05 0.00 0.50 0.00 0.00 0.00 0.00 0.30 0.00 0.00 0.00 
0, U8, 0.30 0.00 0.00 0.00 0.00 0.00 0.20 0.00 0.80 0.80 
03 0.70 0.00 0.00 0.00 0.50 0.00 0.00 0.00 0.00 0.00 


0, U8; 0.00 0.00 0.00 0.00 0.00 020 0.00 0.00 0.00 0.00 
0, U 03 0.00 0.00 0.00 0.20 0.50 020 0.00 0.00 0.00 0.00 
0,U6,U6, 0.00 0.00 1.00 080 0.00 060 000 1.00 0.20 0.20 





4 Numerical Comparison and Discussion 


In this section, we apply the proposed method and the fuzzy DEMATEL 
method to conduct a numerical simulation case, in which the key factors will be 
selected by two methods respectively. Afterwards, the results are compared and 
discussed. 


4.1 The Proposed Method 


Step 1 Define DEMATEL problem and soft CRP parameters 

Similar to the examples in the literatures about DSmT (Smarandache et al. 
2011), here the group IDR matrix construction of the given example works on 
the classical power set 2°, not on the hyper-power set D®. In the numerical simu- 
lation case, suppose the set of factors is F = {f1.f5,f5.f4.fs}, the set of experts is 
E = {é,€,€3,€4,€5}, the set of grade levels is © = {6,,0,,0,} = {0,1,2}, the 
threshold value is 7=0.35, and the acceptable CL threshold values at expert level is 
€,=0.55, and at pair-factors level is €-=0.50, the round counter is t, and the maxi- 
mum number of rounds is ty,y = 5. 

Step 2 Expert assessments extraction and fusion 

Experts are advised to indicate the influence degree to which he or she believes 
factor f; has an effect on factor f; (denoted by f; > f;) under the framework of DSmT. 
As an example, the BBA functions on each pair of factors f; > f; from expert e; are 
shown in Table 1. Due to limited space, other experts’ assessments are not given in 
this paper. 

On the basis of the experts’ assessments, we calculate the similarity matrix by the 
Euclidean similarity function as given in Definition 8 and shown in Eq. (19). 


1.00 0.20 0.16 0.13 0.18 
0.20 1.00 0.15 0.17 0.19 
S=/|0.16 015 100 0.14 0.18 (19) 
0.13 0.17 0.14 1.00 0.18 
0.18 0.19 0.18 0.18 1.00 
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Table 2 Fusion results of expert assessments in the first round 


Ah heh Ah irk heh hoh bok ho hoh >h 


0, 0.00 0.23 0.00 0.13 0.00 0.05 0.46 0.00 0.07 0.05 
05 0.00 0.22 0.39 0.26 0.35 0.20 0.17 0.15 0.25 0.05 
0, UA, 0.08 0.16 0.00 0.00 0.03 0.05 0.08 0.00 0.30 0.32 
03 0.66 0.23 0.19 0.28 0.20 0.58 0.12 0.56 0.22 0.43 


0, UA; 0.02 0.06 0.00 0.00 0.00 0.05 0.16 0.06 0.01 0.06 
0, U 03 0.02 0.00 0.17 0.07 0.20 0.04 0.00 0.05 0.07 0.04 
0, U0, UO; 0.22 0.13 0.26 0.28 0.23 0.06 0.03 0.18 0.09 0.06 





Accordingly, Sup(e,)'! = 0.20 + 0.16 + 0.13 + 0.18 = 0.68,  Sup(e,)! = 0.72, 
Sup(e3)! = 0.64, Sup(e,)! = 0.62, and Sup(es)! = 0.73; the expert weight is com- 
puted as wh! = Sup(e,)'/ X} Suple)! = 0.68/3.40 = 0.20, w! = 0.21, 
w>! = 0.19, w+! = 0.18, and w°! = 0.22; the discounting parameters are derived by 
wh! = wb! /max{w*!|k = 1,...,5} =0.20/0.22 = 0.93, w! = 0.98, w>! = 0.87 
,w*! = 0.82, and w°! = 1.00. Taking be and w! into Shafer’s discounting method 


(as in Eq. (11)), we derive the discounted assessments mi ij=l,...,5, and 


k = 1,2,3,4,5. With the help of the MATLAB code for PCRS, we fuse the BBA 
functions and derive a holism result for the influence degrees between any two fac- 
tors, as shown in Table 2. 

Step 3 Soft consensus reaching process 


k 
The moderator calculates the consensus measure values c 


ï 


be and m; to identify whether an acceptable CL has been reached. Taking expert 


"and ck! by means of 


and fusion assessment BBA functions into Eq. (13) and Eq. (14), we calculate the 
consensus measure values at two levels, as shown in Tables 3 and 4. 


Table 3 The consensus measure values at pair-factors level in the first round 
ich firth Aha fick hoh hoh ho fof foh fo 


ey 0.78 0.65 0.38 0.53 0.58 0.40 0.80 0.29 0.56 0.53 
ez 0.74 0.61 0.50 0.53 0.35 0.42 0.78 0.79 0.56 0.79 
ez 0.78 0.68 0.32 0.57 0.57 0.20 0.56 0.45 0.52 0.55 
e, 0.27 0.57 0.53 0.68 0.37 0.35 0.80 0.38 0.66 0.70 
es 0.78 0.67 0.44 0.57 0.49 0.67 0.57 0.77 0.66 0.60 





The underlined data indicates the pair of factors that do not reach the set consensus level 


Table 4 The consensus measure 
values at expert level in the first 


round chil 0.55 0.61 0.52 0.53 0.62 


e ez ez e4 es 


The underlined data indicates the experts that do not reach the set 
consensus level 
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Table 5 Adjusted assessments of expert e, in the second round 


Ah heh Ah heh heh hoh bok ho hoh >h 


0, 0.00 0.00 0.00 0.50 0.00 0.04 0.00 0.00 0.00 0.00 
05 0.00 0.00 0.22 0.50 0.80 0.16 0.00 0.28 0.00 0.00 
0, U8, 0.30 0.30 0.00 0.00 0.20 0.04 0.20 0.00 0.00 0.70 
03 0.70 0.30 0.37 0.00 0.00 0.46 0.40 0.28 0.00 0.00 


0, U8; 0.00 0.00 0.00 0.00 0.00 0.04 0.40 0.07 0.20 0.30 
0, U 83 0.00 0.00 0.07 0.00 0.00 0.03 0.00 0.07 0.50 0.00 
0, U0, U8; 0.00 0.40 0.34 0.00 0.00 0.24 0.00 0.30 0.30 0.00 





Table 6 Adjusted assessments of expert e, in the second round 


hh Ah Ah heh hoh hoh bok ho h>h >h 


0, 0.00 0.00 0.00 0.00 0.00 0.03 0.50 0.00 0.00 0.00 
05 0.00 0.00 0.00 0.50 0.22 0.14 0.30 0.08 0.00 0.00 
0, U0, 0.06 0.30 0.00 0.00 0.02 0.02 0.20 0.00 0.40 0.00 
03 0.48 0.70 0.50 0.50 0.13 0.39 0.00 0.36 0.60 0.60 


0, U8; 0.02 0.00 0.00 0.00 0.00 0.05 0.00 0.14 0.00 0.20 
0, U 03 0.02 0.00 0.50 0.00 0.13 0.12 0.00 0.12 0.00 0.20 
0, U8, UO, 0.43 0.00 0.00 0.00 0.51 0.27 0.00 0.30 0.00 0.00 


According to the selected threshold value £, = 0.55 and Eq. (16), we find that 
the acceptable CL is not reached in this round for F, # Ø. Then, the moderator lets 
t= 1+1=2, which is lower than fy,,, and carries out the feedback mechanism. 
The order from expert level to pair-factors level should be followed to determine 
the inconsistent experts and the pairs of factors that need to be modified (for those 
values less than £, = 0.50). Tables 3 and 4 show that expert e, greatly differs from 
the group results on the pairs of factors f, > fi, fp > fy, and f} > fs, and e, greatly 
differs from the group results on the pairs of factors fi > fh, fp > fy, fo > fy, and 
J, > fs (see those underlined data in Tables 3 and 4). Experts e, and e, need to 
modify those assessments on the mentioned pairs of factors. Suppose the adjusted 
assessments for expert e, and e, are shown in Tables 5 and 6. 

Because some of the expert assessments are changed. In the second round, expert 
weights are recalculated to be w!? = 0.19, w? = 0.20, w?? = 0.20, w*? = 0.21, and 
w>? = 0.21, and the corresponding discounting parameters are standardized as 
w!2 = 0.89, W? = 0.93, 72 = 0.96, w+? = 1.00, and >? = 1.00. Taking bie and 
ie (Vi, j) and fuse these 
results by PCR5 to construct the group IDR matrix in the second round. The moder- 
ator calculates the new consensus measure values at two levels. Fusion results and 
consensus measures at expert level in the second round are shown in Tables 7 and 8. 


w*? into Eq. (11), we obtain the discounted assessments m 
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Table 7 Fusion results of experts’ assessments in the second round 


hoh hrh Ah irk heh hoh bok ho bhoh >h 








0, 0.00 0.21 0.00 0.12 0.00 0.02 0.44 0.00 0.06 0.06 

05 0.00 0.20 0.48 0.26 0.42 0.10 0.17 0.06 0.22 0.05 
0,U6, 0.16 0.17 0.00 0.00 0.05 0.12 0.08 0.00 0.26 0.29 

03 0.69 0.27 0.13 0.23 0.20 0.53 0.14 0.74 0.28 0.43 

0, U 03 0.04 0.06 0.00 0.00 0.00 0.05 0.17 0.04 0.02 0.07 

0, UO; 0.04 0.00 0.16 0.06 0.27 0.08 0.00 0.02 0.08 0.05 
0,U6,U6, 0.07 O11 0.24 033 0.08 014 O01 015 0.09 0.06 
Table 8 The consensus measure 4 m F m z 
values at the expert level in the 

second round ck? 0.55 0.63 0.66 0.68 0.62 





As shown in Table 8, the acceptable CL has been reached in the second round, which 
means that the current group IDR matrix is a consensus result and can be applied to the 
group DEMATEL method. Then we derive pignistic probability transformation for each 
element in final BBA-formed group IDR matrix M? = [m] sxs: Take the pair of factors 


fı > fa as an example, m,,(0,) = m(0,) + +, U 0,) + (0, U 03) + MO, UO, U 03) = 0.13, 


m2(02) = 0.13, m2(03) = 0.76, and į = 0 x 0.13 + 1 x 0.13 +2 x 0.76 = 1.64. 
After pignistic probability transformation, the final group IDR matrix is constructed as 


in Eq. (20). 


cy 
ll 


0.00 
0.00 
0.00 
0.00 
0.00 


1.64 
0.00 
0.00 
0.00 
1.12 


0.99 
1.33 
0.00 
0.00 
1.25 


Step 4 Calculate the normalized IDR matrix 


Let 


p= 0.00001 


oH Ly 
8, =max(max, <i<p ins Bi M+ MAX] 5; 
matrix can be calculated by Eq. (1). The result is shown as in Eq. (21). 


0.00 
0.00 
0.00 
0.00 
0.00 


0.33 
0.00 
0.00 
0.00 
0.22 


1.22 
1.53 
0.67 
0.00 
0.00 


and 


0.20 
0.27 
0.00 
0.00 
0.25 


0.24 
0.31 
0.13 
0.00 
0.00 


1.13 
0.00 
0.00 (20) 
1.76 
0.00 
we have 


<L a 8) © 4.98, and the normalized IDR 


0.23 
0.00 
0.00 
0.35 
0.00 


(21) 
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Table 9 The attribute 











parameters of factors í É ne — 
fi 0.41 0.00 0.41 0.41 
h 0.19 0.19 0.38 0.00 
Í 0.02 0.24 0.26 —0.22 
fa 0.13 0.23 0.36 —0.10 
Ís 0.13 0.22 0.35 — 0.09 
05ra Coo 
i = factor 
i Jis 
0.25 
, | | | 2 | | rte 
0.1 0.2 0.3 0.4 0.5 
Pay 
"f4 
-0.25 | Js 
-0.5 = 





Fig.3 The cause-effect relation diagram derived by the proposed method 


Step 5 Compute the total relation matrix 
The total relation matrix can be calculated by Eq. (2). The result is shown as 


0.00 
0.00 
A =| 0.00 
0.00 
0.00 


0.14 
0.00 
0.00 
0.00 
0.05 


0.08 
0.08 
0.00 
0.00 
0.08 


0.10 
0.11 
0.02 
0.00 
0.00 


0.09 
0.00 
0.00 (22) 
0.13 
0.00 


According to A = [a;;]sxs and Eq. (3), we derive the following parameters: the 
total influence r; = È a; given by the factor f; to other factors, the total influence 
C= 2a aj; received by the factor f; from other factors, the degree of the impor- 
tant role r; + c; and the net influence r; — c;. These parameters are listed in Table 9. 
Additionally, we construct the cause-effect relation diagram of factors with the hori- 
zontal axis r + c and the vertical axis r — c, as shown in Fig. 3. According to Defini- 
tion 4 and Fig. 3, f, is known as net causer, whereas f3, f,,f; are net receivers. 

Step 6 Obtain major factors with the threshold value 

Because of the threshold value 7 = 0.35, we select the factors whose degrees of 
the important role are greater than the threshold (r; + c; > 4) as major factors. As a 
result, the major factors in the complex system are F = iJa Js} 
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Table 10 The fuzzy linguistic 


en Linguistic terms Triangular fuzzy numbers 
High influence (H) (0.75, 1.00, 1.00) 
Low influence (L) (0.25, 0.50, 0.75) 
No influence (No) (0.00, 0.00, 0.25) 





We find that the Algorithm 1 can be programmed easily. Thus, the proposed 
method in this study is valid and applicable for solving group DEMATEL problems. 


4.2 The Fuzzy DEMATEL Method 


The fuzzy DEMATEL method proposed by Wu and Lee (2007) is a vital and sig- 
nificant improvement of DEMATEL. Since the fuzzy DEMATEL method is not 
only a major extension of DEMATEL but also extracts information with fuzzy lin- 
guistic terms, here we employ it to make a comparison with the proposed method. 
To ensure the results of two kinds of methods can be compared with each other, 
the set of factors F = {f.5,.f3.fa,f5}, the set of experts E = {e}, e2, e3, €4, €s } and 
the threshold value is 7=0.35 are all the same as in Sect. 4.1. Besides, the inputs 
of fuzzy DEMATEL method should be generated from the initial assessments 
given by experts in the proposed method to ensure its comparability. According 
to the procedure of fuzzy DEMATEL method, its computation processes can be 
summarized as follows. 

Step 1 Define DEMATEL problem and fuzzy linguistic scale 

The DEMATEL problem is defined as the same as in Sect. 4.1 and the fuzzy 
linguistic scale is defined as in Table 10, where {No, L, H} are equal to {6,, 02, 63} 
as defined in this paper. 

Step 2 Exact and fuse expert assessments 

In order to ensure the comparability, the inputs of fuzzy DEMATEL method 
should be generated from the initial assessments given by experts in the proposed 
method. Because the assessments given by experts in Sect. 4.1 are in the form 
of BBA functions, how to make a transformation from the BBA functions to the 
inputs of fuzzy DEMATEL method is quite important. It is reasonable and logi- 
cal to choose the grade level with the highest probability as expert assessments 
in the fuzzy DEMATEL method. Following this thought, the inputs of fuzzy 





Table 11 The transformed 


assessments of expert e; fi h fs fa fs 
fi No H L No No 
h No No H H No 
f No No No L No 
fa No No No No H 
ts No No L No No 
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DEMATEL method are transformed from the experts’ assessments in Sect. 4.1. 
For example, the transformed assessments of expert e, are shown in Table 11. 

As introduced in fuzzy DEMATEL, the Converting Fuzzy data into Crisp 
Scores (CFCS) defuzzification method is applied to aggregate these assessments 
by five experts. Due to limited space, the detailed CFCS steps can be referred to 
Wu and Lee (2007) and are not repeated here. The IDR matrix G’ = [Silsxs is pro- 


duced as in Eq. (23). 


0.04 0.19 0.87 0.19 0.50 
0.19 0.04 0.96 0.96 0.04 
G' =| 0.04 0.04 0.04 041 0.19 (23) 
0.04 0.04 0.04 0.04 0.87 
0.04 0.50 0.59 0.04 0.04 


Step 3 Calculate the normalized IDR matrix 


The normalizing method for IDR matrix in fuzzy DEMATEL is the same as 
in traditional DEMATEL. The normalized IDR matrix is calculated by taking 
Eq. (23) into Eq. (1) and it is shown as in Eq. (24). 


0.02 0.08 0.35 0.08 0.20 
0.08 0.02 0.38 0.38 0.02 
D' =| 0.02 0.02 0.02 0.16 0.08 (24) 
0.02 0.02 0.02 0.02 0.35 
0.02 0.20 0.24 0.02 0.02 


Step 4 Compute the total relation matrix 
The computing method for total-relation matrix in fuzzy DEMATEL is also the 
same as in traditional DEMATEL. The total-relation IDR matrix is calculated by 


taking Eq. (24) into Eq. (2) and it is shown as in Eq. (25). 





0.02 0.01 0.18 0.02 0.07 
0.01 0.02 0.21 0.20 0.00 
A’=|0.00 0.00 0.02 0.03 0.01 (25) 
0.00 0.00 0.00 0.02 0.14 
0.00 0.05 0.09 0.00 0.02 
Table 12 The factors’ attributes F E ree Fed 
parameters 
fi 0.30 0.03 0.33 0.27 
h 0.44 0.08 0.52 0.36 
f 0.07 0.50 0.57 — 0.43 
Ja 0.17 0.28 0.44 —0.11 
Ís 0.16 0.25 0.41 —0.09 
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Fig.4 The cause-effect relation diagram derived by fuzzy DEMATEL method 


According to Eq. (3), we derive the following parameters as shown in Table 12 
and the cause-effect relation diagram as shown in Fig. 4. According to Definition 
4 and Fig. 4, it is obvious to find that fi, fọ are net causers, while f3, f4, fs are net 
receivers. 

Step 5 Set a threshold value and obtain the major factors 

Because of the threshold value 7 = 0.35, only the factors with r; + c! > 7 should 
be chosen as the major factors. As a result, the major factors in the complex system 


are Ê! = {f,, fa, fasts} 
4.3 Discussion 


It is obvious to find that the major factors in the complex system determined by the 
proposed method are Ê = {fi -So-fa.f5}, while those determined by the fuzzy DEM- 
ATEL method are F’ = {f;,f;,f,,f5}- The results of two kinds of methods are differ- 
ent from each other. Which one is more reasonable? Now we make discussions from 
the following three aspects. 

(1) Expert assessment extraction mechanism. The assessments of experts are 
used as the fundamental inputs of DEMATEL whether in the proposed method 
or in the fuzzy DEMATEL method. The proposed method allows experts to give 
assessments with BBA functions, while the fuzzy DEMATEL method employs 
the fuzzy linguistic scale to express their assessments. The local or global igno- 
rance in experts’ minds can be well reflected in the proposed method (e.g., 
be = {((0), 02), 0.30), (03, 0.70) }), but it is unfortunate to find that such the igno- 
rance cannot included in the fuzzy DEMATEL method (it seems that only the 
grade level/fuzzy linguistic scale with the highest probability may be allowed to 
express assessments, e.g., bi = {L} = {6,}). Consequently, it is believed that 
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the proposed method is more feasible than the fuzzy DEMATEL method in the 
aspect of expert assessment extraction mechanism. 

(2) The importance roles of experts. As discussed in Sect. 3.1, the importance 
roles of experts are reflected by expert weights that are calculated based on 
expert assessments information, hence, the weight parameters could effectively 
express the relative importance of experts in the group and they have a signifi- 
cant impact on the group DEMATEL decision results. If the importance roles 
of experts are neglected in the process of decision making, the decision results 
may lose effectiveness. The proposed method calculates expert weights based 
on expert assessments with the aid of evidence distance and employs Shafer’s 
discounting method to modify the subjective assessments. Unfortunately, the 
importance roles of experts are unconsidered in the whole processes of fuzzy 
DEMATEL method. Consequently, it is believed that the proposed method is 
more accurate than the fuzzy DEMATEL method in the aspect of reflecting the 
importance roles of experts. 

(3) Group consensus reaching. Group consensus reaching is a key problem in 
the GDM field. Only when the assessments given by experts reach an acceptable 
CL, the GDM results are seen to be valuable and credible. In other worlds, the 
GDM results that lack of consensus may be ineffective and have few reference 
value for decision-making. From Sect. 4.2, we know that the fuzzy DEMATEL 
method is a static method without taking the group consensus reaching into con- 
sideration. In the proposed method, we apply the soft CRP into the construction 
of group IDR matrix to help the experts group reach a high CL, among which 
the feedback and modification mechanism are introduced. Consequently, it is 
believed that the proposed method is more reasonable than the fuzzy DEMATEL 
in the aspect of reaching group consensus. 


5 Conclusions 


In the present study, the DSmT is used to extract and fuse expert assessments, 
and a soft CRP is introduced to construct the group IDR matrix for the DEMA- 
TEL. The DSmT-based group DEMATEL method and the corresponding algo- 
rithm are proposed. Moreover, a numerical comparison is performed to discuss 
the applicability of the proposed method and algorithm. The main contributions 
of the present study can be summarized into three aspects. 

Firstly, an expert assessment extraction and fusion mechanism is established 
on the basis of DSmT. Expert assessments on the influence relations among fac- 
tors are extracted by BBA functions, which can help experts to express uncer- 
tainty and incompleteness assessments. The PCR5 of DSmT, which can overcome 
the defects of intuitional paradox in Dempster’s combination rule, is employed to 
make fusion for the individual BBA functions discounted by Shafter’s discount- 
ing method. 

Secondly, expert weights are defined and introduced to reflect importance roles 
of experts in a group. Following the principle of pairwise comparisons, expert 
weights are calculated by the Euclidean similarity function to reflect experts’ 
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relative importance in the group. Expert weights are not fixed during the whole 
processes of group IDR matrix construction and they may vary dynamically with 
different expert assessments in each round. Shafer’s discounting method is used 
to discount expert assessments in each round so as to reflect the importance roles 
of experts in the group dynamically. The above processes are beneficial to obtain 
an accurate group IDR matrix for the DEMATEL. 

Thirdly, a soft CRP is established to construct the group IDR matrix with 
consensus and an algorithm is summarized for the group DSmT-based DEMA- 
TEL. The consensus measures at expert and pair-factors levels are defined to help 
establish the feedback/modification mechanism, based on which a soft CRP is 
established for the construction of consensus group IDR matrix, so that experts 
can reach a high CL. An algorithm for DSmT-based group DEMATEL method 
with reaching consensus is proposed to identify major factors in a complex sys- 
tem. The proposed algorithm can be programmed easily and is valid and applica- 
ble for solving group DEMATEL problems. 

The proposed method forms an expert assessment extraction mechanism on the 
basis of the BBA function. However, the belief degrees or probabilities in BBA 
functions given by experts may hardly be assessed with exact values in more 
complex situations. Therefore, investigating how to deal with the group DEMA- 
TEL with the interval BBA function may be a good direction for future research. 
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